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The excitons in the orthorhombic phase of the perovskite CHsNHsPbls are studied using the 
effective mass approximation. The electron-hole interaction is screened by a distance-dependent 
dielectric function, as described by the Haken potential or the Pollmann-Biittner potential. The 
energy spectrum and the eigenfunctions are calculated for both cases. The effective masses, the low 
and high frequency dielectric constants, and the interband absorption matrix elements, are obtained 
from generalized density functional theory calculations. The results show that the Pollmann-Biittner 
model provides better agreement with the experimental results. The discrete part of the exciton 
spectrum is composed of a fundamental state with a binding energy of 24 meV, and higher states that 
are within 2 meV from the onset the unbound exciton continuum. Light absorption is dominated by 
the fundamental line with an oscillator strength of 0.013, followed by the exciton continuum. The 
calculations have been performed without fitting any parameter from experiments and are in close 
agreement with recent experimental results. 
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Hybrid organic/inorganic perovskites based on lead or 
tin tri-halides are semiconductor materials that have rev¬ 
olutionized the research of thin film solar cells. With the 
first prototypes demonstrated six years agOji record cell 
efficiencies have surpassed the barrier of 20%i^‘^ Methyl- 
ammonium lead iodide (CHaNHsPbla) is one of the most 
studied members of this family, and it has been applied 
as photon absorber and charge transporting material)^ 

CHaNHsPbla presents two phase transitions at ^ 
162 K and ^ 327 K. At these transitions, the crystal 
symmetry changes first from orthorhombic to tetrago¬ 
nal and then to cubic symmetryj^^— The three phases 
differ by small changes of the lattice vectors, rotations 
of the characteristic Pble octahedra, and the orienta¬ 
tion of the CHsNH;^ cations. In the tetragonal and cu¬ 
bic phases, the CHsNH;^ cations present orientational 
and dynamic disorder,— with a deep effect on the dielec¬ 
tric properties— In the low temperature orthorhom¬ 
bic phase, the CHaNH;^ positions and orientations are 
fixedidiid 

The electronic band structure of CHaNHaPbla has 
been explained on the basis of generalized density func¬ 
tional theory (hybrid functionals) or Green functions 
GW calculations, in both cases including the spin-orbit 
coupling— yid For the orthorhombic phase, the valence 
band maximum (VBM) and the conduction band min¬ 
imum (GBM) are located at the T point corresponding 
to the 48-atoms unit cell, and the fundamental gap is 
1.68 eV— Both the VBM and GBM are doubly degen¬ 
erated, with nearly symmetric effective mass tensors. 

Exciton peaks are observed in the light absorption 
spectra at low temperaturerdr— just below the inter¬ 
band absorption edge, or melded with it, depending 
on the temperature. According to the Wannier-Mott 
modelf^i^i the exciton is similar to a hydrogen atom with 
the proton and electron masses replaced by the hole and 


electron effective masses, and the Coulomb interaction is 
screened by a dielectric constant e. Therefore, the exciton 
binding energy and the Bohr radius are Ry = /ie^/2fi^e^ 
and aex = h^e/y.e^, where yt = merrih/(me + mh) is the 
reduced electron-hole mass. 

One distinct feature of CHsNHaPbIs is the large dif¬ 
ference between the static dielectric constant cq and the 
high frequency constant Coo, i-e., for frequencies higher 
than those of the phonon absorption. Values of Coo in 
the range 4.5 — 6.5 have been calculatedwhile 
values close to 25 have been estimated for eg— Such 
difference is larger than in traditional inorganic semi¬ 
conductor and should cause important polaron effects, 
such as the effective mass and gap renormalization, as 
well as and non-hydrogenic exciton states. For the lat¬ 
ter, immediately arises the question wether the screening 
constant e should be the static dielectric constant eg or 
the high frequency Coo- Using the values listed in Table 
m the static and the high frequency dielectric constants 
lead to very different values of the exciton binding en¬ 
ergy Ryo = 2.8 meV and Ryoc = 50 meV, respectively. 
Such different energies lead to different conclusions with 
respect to exciton dissociation due to thermal excitation, 
as well as to different interpretation of luminescence and 
transport properties. 

Early estimations of the exciton binding energyiid^ 
^ 37 — 50 meV were based on measurement of the exciton 
diamagnetic coefficient and interpretation based on the 
hydrogenic model with screening by Coo = 6.5. Recent 
studies of the temperature dependence of photolumines¬ 
cence spectra^l and numerical analysis of the absorp¬ 
tion spectra^ii^ have provided updated exciton binding 
energies around 16-19 meV. The latter values point to 
a screening constant intermediate between Coo and eg. 
Even et al2i htted the absorption spectrum using the 
Wannier-Mott exciton model and obtained an effective 
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dielectric constant eeff = 11. 

In fact, the differences between eg and Coo express the 
electric polarization associated to the optical phonons 
and the electron-phonon interaction. The stationary 
states are coupled states of electronic and the vibra¬ 
tional phonon field. The quantum calculation of these 
coupled states is beyond the current capabilities of ab 
initio methods. Model Hamiltonians^Sr— allow one to 
map the coupled electron-phonon excitations into effec¬ 
tive electronic states, and to obtain the energies of sta¬ 
tionary states. Even when simplifying approximations 
are inherent in the models, they can provide a criterium 
on the relevant dielectric screening constants. In this Ar¬ 
ticle, we apply the model Hamiltonians of Haken^SiSl and 
that of Pollmann and Biittner— to the exciton spectrum. 
This formalism is applicable to the low temperature or¬ 
thorhombic phase because in the tetragonal phase the 
static dielectric increases strongly, associated to the re¬ 
orientation of CHgNH;^ cations, and the exciton effects 
practically disappear 

The strength of the interaction of electrons and optical 
phonons is given by the coupling constant 

ap = ^me^/2h?elELo, ( 1 ) 

where E^o is the energy of the longitudinal optical 
phonon. This model was developed for simple crys¬ 
tals that display one single LO phonon branch. For 
this application, we have chosen Elo as the shift of 
the main peak in the CHaNHaPbIg Raman spectrumi^ 
The ionic screening parameter appearing in Eq. o is 
1/e^ = l/coo - 1/eo- 

For transport properties, relevant after exciton dissoci¬ 
ation, polaron masses must be considered rather than the 
bare electronic masses computed with fixed ions. They 
can be estimated using the Frohlich’s continuum theory 
of the large polaron^^S which predicts 

m = TO (^1 -I- j . 

The polaron bands undergo an additional shift given by 
AEp = —apE^o- With the data of Table HI this leads to 
a reduction of the electronic band gap by 95 meV. 

The Haken mode P^i^^ describes two interacting po- 
larons, each one with a radius much smaller than the 
exciton effective radius, and expresses the effective po¬ 
tential for the electron-hole Coulomb interaction as 



Here le,h = \/W/2m^^hMio are the electron- and hole- 
polaron radii determined using bare band electron and 
hole effective masses. The polaron effective mass param¬ 
eters must be used in the kinetic energy terms of the 
Hamiltonian 

The model proposed by Pollmann and Biittner— (PB) 
takes into account the correlation between electron and 



FIG. 1. Haken and PB potentials compared with the Coulomb 
potential screened by eo and too. Also shown are the eigenen- 
ergies for each potential. 


hole polarons, and leads to corrections to the Haken po¬ 
tential. The resulting electron-hole interaction potential 
is 


VpB{r) 


egr e*r VAto Ato / ’ 


with Ato = TO/j — TOg. This potential was derived assum¬ 
ing le,h that the polaron lengths are much smaller than 
the effective exciton radius, which entered as a variational 
parameter in the original calculations.— The hare band 
electron and hole masses must be used in the kinetic en¬ 
ergy terms of the PB Hamiltonian. 

In the present work, the exciton energies are obtained 
solving the radial Schrodinger equation for the relative 
coordinate wave function of the exciton $(r) 


d2$ 

dr^ 


2d$ /2u , , + 


where V{r) is the electron-hole interaction potential 
(Coulomb, Haken, or PB), I is the azimuthal quantum 
number, of which we only consider I = 0 that are the opti¬ 
cally active states. The Eq. (g]) for Z = 0 has been solved 
integrating the equation starting from r = 0 with the 
conditions $(0) > 0, $'(0) = 0 and imposing $(rc) = 0, 
where rc is a cutoff radius sufficiently large to mimic the 
boundary conditions at infinity. The cutoff radii Vc are es¬ 
tablished solving the equation for the Coulomb potentials 
and comparing the numerical energies with the known ex¬ 
act solutions. We have used exciton atomic units og and 
Ryo for the radius and energy, respectively. The func¬ 
tions are normalized according to 

c 

47r|$(r)|^r^dr = 1. (5) 

The optical oscillator strengths are defined as 

/n = |^|(^n|C-^10)P, (6) 

^n,0 
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TABLE I. Parameters defining the polarons in CHsNHsPbls. 


Dielectric constants 


Coo 

5.32 

Coo/eo 

0.236^ 

LO phonon energy 


Elo 

38.5 meV^ 

Coupling constants 


ae 

1.18 

Cth 

1.28 

Bare carrier masses 


melmq^ 

0.190 


0.225 

Polaron masses 


ml/mo 

0.228 

ml/mo 

0.273 

Polaron radii 


le 

22.83 A 

Ih 

21.00 A 

Polaron shift 


ae; 

-45.3 meV^ 

ae; 

-49.2 meV^ 



where mo is the free electron mass, huJn,o = Eg* + is 
the transition energy, |0) and |iE'„) are the ground and ex¬ 
cited states of the crystal, respectively, and v = i[H, r\/h 
is the velocity operatori^ E* is the renormalized gap 
(with the polaron shift), and En are the eigenvalues of 
Eq. Q. We shall approximate /„ by the expression for 
pure excitons, i.e., neglecting the phonon coupling, as 


Therefore, we obtain the simplified expression 


fn 


^Ucv 

hUnfi aj^ 




( 9 ) 


Let us stress that the exciton Bohr radius Uq appears in 
Eq. ([9]) only if the normalization condition (O is applied 
in relative units of Oq. 

Both the Haken and Pollmann-Biittner potentials be¬ 
have like a Coulomb potential for very large distance 
{r ^ lejh) or very short distances (r <C lejh), screened 
by the low and high frequencies dielectric constants, re¬ 
spectively. Figure [T] shows in logarithmic scale, the limit¬ 
ing Coulomb potentials screened by ep and Coo- These are 
represented by the straight lines, enclosing the Haken and 
Pollmann-Biittner potentials, that interpolate the limit¬ 
ing cases. Horizontal lines represent the eigenenergies 
of the exciton relative motion. The axes in the figure 
are in units of static (fully screened) exciton radius and 
Op = h^eo/ge^ and exciton energy Ryo = In 

these units, the static Coulomb potential is given by —2/r 
and the exciton eigenenergies are E^ = —Xjr?. The 
Coulomb potential and the hydrogenic energies defined 
by Coo are where = ep/coo- For the 

parameters of CHaNHaPblp (cr = 4.25), one can appre¬ 
ciate in Figure [1] and Table [III that the lowest exciton 
levels are E^ = —11.26 and E^^ = —8.65 for the Haken 
and PB potentials. These values represent a significant 
correction to either E^ = — 1 or E^ = —18. The excited 
exciton energies of Haken and PB potentials approach 
the values —1/n^ for high n. 

In order to compare the energies E^ and E^^ one 
must consider that Ryo is defined either by the po¬ 
laron or the bare reduced mass in the first and second 
model, respectively. In absolute units, E^ = —37 meV 
and EC^ = —24 meV. Is seems that the PB value is 
in better agreement with the experimental values near 
19 meV,^i25 


fn 


2mp 

fujn,0 


E E ^IKo|i)c|«.o)P%^|4>„(0)r 

O QjQ 

cv a—x,y,z ^ 

( 7 ) 


In the above expression, H/.u, is the normalization vol¬ 
ume of the center-of-mass part of the exciton envelope 
wave function function, which we consider as the volume 
of one formula unit, i.e., one fourth of the unit cell volume 
952.5 A^. With this convention, the oscillator strength 
is equivalent to the values reported elsewhereJ^d^ The 
factor 1/3 and the sum in a correspond to isotropic av¬ 
erage of the crystal orientations, u^q and Uco are the 
Bloch functions of the valence band maximum and con¬ 
duction band minimum, which in this case are both dou¬ 
bly degenerate. Using first principles calculations (see 
the Appendix) we have calculated the parameter 


100 


I 10 

I& 

m 

a 


1 

0 2 4 6 8 10 12 14 16 

n (level) 

FIG. 2. Oscillator strengths of polaron excitons and hydro¬ 
genic exciton. 
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C7cn = ^E E ^l('aco|i)a|Mr;o)P = 1-706 eV. (8) 

cv a—x,y,z 


The values of $(0) in Table |TT] and Figure [5] shows 
that the ground exciton is optically active. According to 
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TABLE II. Exciton binding energies and wave function at 
origin. The wave function is normalized with the radii in 
units of oq. In the same units, the hydrogenic functions fulfill 
4?n(0) = l/VTrn^. In normal units, 4>„(0) must be divided by 


n 

Haken model 
E^/Ryo (0) 

PB model 

E^^/Ryo K^iO) 

1 

-11.2605 

4.7255 

-8.6473 

4.5421 

2 

-0.8128 

1.0832 

-0.4622 

0.6507 

3 

-0.2092 

0.3636 

-0.1661 

0.3031 

4 

-0.0979 

0.2047 

-0.0839 

0.1811 

5 

-0.0567 

0.1358 

-0.0505 

0.1236 

6 

-0.0370 

0.0985 

-0.0337 

0.0911 

7 

-0.0260 

0.0756 

-0.0240 

0.0708 

8 

-0.0193 

0.0604 

-0.0180 

0.0570 

9 

-0.0149 

0.0497 

-0.0140 

0.0472 

10 

-0.0118 

0.0418 

-0.0112 

0.0399 

11 

-0.0096 

0.0358 

-0.0092 

0.0343 

12 

-0.0080 

0.0311 

-0.0076 

0.0299 

13 

-0.0067 

0.0274 

-0.0065 

0.0264 

14 

-0.0057 

0.0243 

-0.0055 

0.0235 

15 

-0.0050 

0.0218 

-0.0048 

0.0211 


Eq. © the oscillator strength in the PB model is 0.013. 
The second exciton level is only 2.7 meV (Haken) or 
1.2 meV (PB) below the edge of the continuum spectrum, 
and their oscillator strengths are one order of magnitude 
smaller than for the main line (n = 1), making these tran¬ 
sitions practically undetectable in the optical spectra. 
Higher energies approach to the sequence Eg — Ryo/v?. 
For these higher levels, the oscillator strengths become 
proportional to the hydrogenic oscillator strengths. The 
proportionality constant is fitted to Ph = 4.282 (Haken) 
and Ppb = 4.214 (PB) in each case. The fitting function 
was /(n) = P + j/{n — no), and the fitted P are stable 
for any subset of data with n > 7. Let us note that 

P^^. ( 10 ) 

€oo 


This result, together with the approximation of the en¬ 
ergies by the sequence E*—Ryo/n^, leads to to a constant 
absorption spectrum near the band gap energy similar to 
the case of hydrogenic exciton^i^ 


^ EgUrCm^ 


( 11 ) 


where is the refraction index and can be approximated 
by y/f-ca- Evaluation the physical constants, and using 
Eqs. (fTOl) . ([8]), and (IA5I) . the above expression can be 
cast as 




(3.466 nm ^) 


C/e. 


Ur Eg 



( 12 ) 


Using the material parameters of CHaNHaPbla, a{Eg) ^ 

3.1 

The exciton spectrum obtained presents of a low en¬ 
ergy non-hydrogenic state that is 37 or 24 meV (in Haken 
or PB models, respectively) below the onset of continuous 
free polaron spectrum. The rest of the discrete exciton 
spectrum resembles the hydrogenic series with static di¬ 
electric constant eo. The first exciton state dominates the 
absorption edge with oscillator strength ^ 0.013, while 
the oscillator strengths of higher states are one or two 
orders of magnitude smaller. Like the case of hydrogenic 
exciton, the coalescence of absorption lines from higher 
states generates a quasi-continuum spectrum with a fi¬ 
nite value at the gap energy. The absorption coefficient 
at gap energy in enhanced with respect to the hydrogenic 
model by a factor eo/coo- Extrapolating this behavior to 
the room temperature tetragonal phase, this enhance¬ 
ment factor could be responsible of the high capacity of 
this material to collect the photon energy in the photo¬ 
voltaic cells. 

To interpret the absorption experiments and to deter¬ 
mine de band gap, one needs to know wether the onset of 
the continuous absorption spectrum corresponds to the 
exciton continuum spectrum, or to the accumulation of 
discrete lines below the band gap. In other words, what is 
the energy range of constant absorption coefficient given 
by Eq. (HU. Considering that the higher exciton levels 
are within 2 meV of the continuum, and that the exciton 
absorption spectrum is dominated by the fundamental 
state, one can conclude that the band gap coincides with 
the absorption threshold after filtering the first exciton 
peak. On the other hand, the polaronic effect downshifts 
the gap by 95 meV. Therefore, the measured gap 1.68 eV 
should be understood as an electronic gap of 1.78 eV de¬ 
creased by the polaron shift. 

Several kinds of estimations of the exciton binding 
energy in CHaNHsPbU have been reported. The first 
method was employed by Hirasawa et alj^ and repeated 
later with improved accuracy by Tanaka et al^ They 
measured the exciton diamagnetic coefficient in magneto¬ 
absorption spectra, and related the measurements with 
the binding energy in the framework of the hydrogenic 
model with the high frequency dielectric constant. With 
this model, Tanaka et al determined a binding energy of 
50 meV. The use of the high frequency dielectric con¬ 
stants was a choice of the model, and not determined by 
the experiments. Let us note that Tanaka et al, and Hi¬ 
rasawa et al used a value Coo = 6.5 that is higher than 
our value. If our value Coo = 5.32 were used, a bind¬ 
ing energy of 65 meV would be obtained. Conversely, 
our binding energies using the larger dielectric constant 
would be smaller. These values of the exciton binding 
energy are strongly biased by the choice of the dielectric 
constants. 

The second method has been applied by Sun et 
{Eb = 19 meV). They obtained the binding energy 
by fitting the photoluminescence intensity as a function 
of temperature with an Arrhenius equation, not using 
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any model of the exciton states. It is rewarding that 
our calculation of the binding energy is in good agree¬ 
ment with the values determined by Sun et al. Huang 
and Lambrecht^^ have argued, in a study of cesium tin 
halide perovskites, that the photoluminescence temper¬ 
ature dependence just give information on the free ex¬ 
citon linewidth or the binding energies of bound exci- 
tons, but not on free excitons. However, Even et al2i 
fitted the absorption spectrum using the Wannier-Mott 
exciton model and obtained a similar value for the bind¬ 
ing energy, and reported an effective dielectric constant 
Eeff = 11. This value of the dielectric constants, to¬ 
gether with the assumed reduced mass O.lGmcri means 
a binding energy of 19 meV. Another method indepen¬ 
dent of the dielectric function has been used by Miyata 
et al^ who performed magneto-absorption experiments 
with very high magnetic fields, determining a value of 
16 ± 2 eV. It is interesting that Miyata et al were able 
to detect the 2s exciton state for high magnetic field and 
extrapolated a Is — 2s difference of 15 meV at low mag¬ 
netic field. Henceforth, assuming the hydrogenic model, 
they estimated the binding energy in 20 meV. However, 
extrapolating the Landau levels of the free exciton spec¬ 
trum they obtained the precise value of 16 meV. This 
observation agrees with our result that the 2 s state is 
within 2 meV of the free exciton edge. 

The evaluation of the oscillator strength provides an¬ 
other argument against the model of hydrogenic Wannier 
excitons screened by cq. Ishiharapi^ reported an experi¬ 
mental values of 0.02, which is close to our value 0.013. 
If the ground exciton were well described by hidrogenic 
model with cq) then <I>i(0) = l/v^ (compare with Ta¬ 
ble H] and Fig. 12), and the oscillator strength would be 
^ 64 times smaller. Therefore, at least for low tempera¬ 
ture, the fundamental exciton state does not correspond 
to screening by Cq. 

We wish to stress that we have not fitted any param¬ 
eter in this work, which would bring the exciton binding 
energies in closer agreement with the recent experimental 
results. The parameters with larger uncertainty are the 
dielectric constants and the LO phonon energy. The only 
available experimental value of Coo = 6.5^ is larger than 
the ab initio value used here and that value would reduce 
the calculated binding energy. The measurement of Coo is 
rather old, with few published details, and a new deter¬ 
mination for present-day thin films would be welcomed. 
The LO phonon energy Elo has been chosen from the 
more prominent peak in the calculated Raman spectrum 
of CHaNHaPblar^ which is close to LO phonon ener¬ 
gies in H-VI and HI-V semiconductors. As mentioned 
above, the model Hamiltonians were developed assuming 
a unique LO phonon energy. The Raman spectrum of 
CHsNHaPbIs shows bands at lower wave numbers. Us¬ 
ing and average energy of the Raman active peaks, which 
do not have necessarily LO character, may lead to lower 
exciton binding energy. An extension of the PB Hamil¬ 
tonian to include several LO phonon branches would be 
a better founded approach. 


In summary, we have calculated the exciton bind¬ 
ing energies and oscillator strengths using two model 
Hamiltonians of the exciton-phonon coupled system. The 
Pollmann-Biittner model Hamiltonian gives a binding en¬ 
ergy in good agreement with recent experimental deter¬ 
minations. The calculated oscillator strength of the main 
exciton line agrees with the value estimated from experi¬ 
ments, while the strengths of higher transitions are much 
smaller. 
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Appendix A: Interband matrix element 

With the VASP code^ the tensor dielectric function is 
computed in the longitudinal approximation^! 


,( 2 ) 


^ lim ^ gswk(5(eck-fq - e«k - 
it q^o 

c,v,]<. 


^ ('^ck+eQ,g I'^Lik) ('^ck+e^g I'^Lik) 


(Al) 


where Wk are the k-point weights, defined such that they 
sum to 1, Cck-i-q are band energies, H is the unit cell vol¬ 
ume, rriQ is the vacuum electron mass, are polarization 
vectors. The factor gs is the spin degeneracy, which is 2 
in Ref. and the bands v, c in the sum are restricted 
to have the same spin. In the calculation with spin-orbit 
coupling, we consider gs = 1 and the sum is over all 
pairs of valence and conduction bands. The Eq. lAll is 
equivalent to the transverse approximation^ 


Sa/3(^) = 2 gSWwSiCcW “ e«k “ Sw) 




c,L;,k 


/ I I \ / I /s TlIvR I , ^ / K n\ 

^ H I'^uk) \^ck|'^/3 H |‘^Dk) ■ (-^^j 

mo mo 


For a local Hamiltonians with spin-orbit coupling, 
mov = —ih\7 + {h/imoc^)<7 x VU. The PAW poten¬ 
tials and the hybrid functionals introduce non-locality in 
the Hamiltonian, and the velocity operator contains addi¬ 
tional termsi^ For the purpose of the optical properties 
of the exciton, we only need the values of (itckl'C/sIw^k) 
between the VBM and the GBM, which occurs at the 
r point {kp = 0). These values (ucklvisluyu) can be fit¬ 
ted from the dielectric function, which is calculated using 
(EH). Hence, if the contribution of the T point can be 
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separated from the other k-point contributions, we have The parameter Pcv has dimension of momentum, and 
that P^^j2mQ is the parameter U^v defined in Eq. ([5]). 


e 


( 2 ) 

ap\T 


(w) 




gswo ^ (5(ec0 - e„o - Sw) 

CV 


X (McO I Wa I UyO ){Uco\v/3\Uyo)* . 


(A3) 


In the above expression, the sum is restricted to the 
top valence bands and bottom valence bands. 

To fit with the exciton, we consider the averaged di¬ 
electric function e = (l/3)Treo,a 


-( 2 ) 


(w) = 




rgswoPcySieco - evo - fihj), (A4) 


with 


1 

3 


/ 

EE \{Uco\moVa\Uyo)\'^. 

OL CV 


(A5) 


In practical calculations, S{eco — (vO — is replaced 
by a smearing function. If gaussian smearing is used 
for the self-consistent calculation, the computed spec¬ 
trum must be fitted with a Gaussian function weighted 
by iTr^e^gswoPcy/^c^^mg. The ab initio calculation was 
performed sampling the Brillouin zone with a 2 x 2 x 2 
k-point grid centered at the T point, which is sufficient to 
obtain total energies and charge densities, but it is coarse 
for calculation of optical properties and it allows to iso¬ 
late the contributions of the T—point transitions. With 
this k-point mesh, the weight wn = 1/8. The details of 
the calculation are given in Ref. Ila 
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